Supercollision cooling in undoped graphene 



A. C. Betz/'Q S. H. Jhang/'B E. Pallecchi, 1 ' 2 R. Feirrera, 1 G. Feve, 1 J.-M. Berroir, 1 and B. Plagais^Q 

'Laboratoire Pierre Aigrain, ENS-CNRS UMR 8551, 
Universites P. et M. Curie and Paris-Diderot, 24, rue Lhomond, 75231 Paris Cedex 05, France 
2 Laboratoire de Photonique et Nanostructures, CNRS-UPR20 CNRS, 
Route de Nozay, 91460 Marcoussis Cedex, France. 



o 

(N 
^— > 
O 

O 

in 



T3 

o 
o 



>: 
ov 

00' 

O' 
(N 



X 



Carrier mobility in solids is generally limited 
by electron-impurity or electron-phonon scatter- 
ing depending on the most frequently occurring 
event. Three body collisions between carriers 
and both phonons and impurities are rare; they 
are denoted supercollisions (SCs) 1 . Elusive in 
electronic transport they should emerge in re- 
laxation processes as they allow for large energy 
transfers 2 . As pointed out in Ref. this is 
the case in undoped graphene where the small 
Fermi surface drastically restricts the allowed 
phonon energy in ordinary collisions. Using elec- 
trical heating and sensitive noise thermometry 
we report on SC-cooling in diffusive monolayer 
graphene. At low carrier density and high phonon 
temperature the Joule power P obeys a P oc T 3 
law as a function of electronic temperature T e . 
It overrules the linear law expected for ordinary 
collisions which has recently been observed in re- 
sistivity measurements 3 . The cubic law is char- 
acteristic of SCs and departs from the T 4 depen- 
dence recently reported for metallic graphene 4 
below the Bloch-Griineisen temperature. These 
supercollisions are important for applications of 
graphene in bolometry and photo-detection. 

Understanding how two-dimensional (2D) electrons in 
graphene relax their energy to the lattice is not only a 
central problem in condensed matter physics but also an 
important issue in the design of graphene devices^. Due 
to the large optical phonon energy of fi ~ 200 meV in 
graphene, the emission of acoustic phonons is the only 
efficient cooling pathway for hot electrons below the en- 
ergy Q. Unlike conventional metals with large Fermi 
surfaces, where the Debye temperature sets the bound- 
ary between high- and low-temperature behavior in the 
electron-phonon interaction, a new characteristic temper- 
ature arises in graphene, the Bloch-Griineisen tempera- 
ture (Tbg)- It results from the small Fermi surface in 
graphene and is defined by the maximum phonon wave 
vector q max —2kp (Fig.QJ}). Above Tbg, only a fraction 
of acoustic phonons with wavevector q <2kp can scatter 
off electrons. While this phase space restriction leads to 
the crossover behavior of electron-phonon resistivity be- 
tween high-temperature p(T) cx T and low-temperature 
p(T) oc T 4 dependence 3 -, it also imposes a significant con- 
straint on the cooling of hot electrons; the energy dissi- 
pated by acoustic phonons cannot exceed As Tbg per one 
scattering event (Fig.QJ:). As a result, the electron- lattice 
cooling can be slow for phonon temperatures T p h 3> Tbg, 



requiring many scattering events to dissipate the heat of 
the hot electrons. 

On the other hand, a recent theory predicts alternative 
cooling pathways mediated by supercollisions 1 . Here, 
disorder-assisted phonon scattering or two-phonon scat- 
tering events enable the emission of acoustic phonons 
with wave vectors q > 2fcp , thus transferring higher en- 
ergy than normal collisions. Such collisions occur less 
frequently than ordinary collisions, but dominate the 
electron-lattice cooling by utilizing the entire thermal dis- 
tribution of phonons (Fig.[TJ;). 

A unique feature of graphene is the tunability of Tbg 
with Fermi energy Tp 3 . As both the electron and the 
phonon energy are linear in wave vector, the Bloch- 
Griineisen temperature is simply given by /cbTbg = 
(2v s /vf)Ef, where v s and vf are the sound and the 
Fermi velocity, respectively. Hence, both energy scales, 
^bTbg and Ep, are linked by the ratio 2v s /vp ~ 0.04. 

In this work, we experimentally test the different mod- 
els of electron-phonon coolingiZr— . The low- (T p h < 
Tbg) and high-temperature (T p h > Tbg) regimes be- 
come accessible by tuning the size of the Fermi surface, 
i.e. tuning Tbg by means of an electrostatic gate poten- 
tial V g (Fig.[7J). Electrons are heated with a bias voltage 
V in a two-probe configuration, and the electron-lattice 
cooling rate is investigated by means of Johnson noise 
thermometry, a primary temperature measurement tech- 
nique based on the fundamental properties of thermal 
fluctuations in conductors. The electron temperature is 
deduced from radio frequency shot noise measurements, 
relying on the relation Si — 4ksT e /R between the cur- 
rent noise spectrum Si, the sample resistance R, and T e . 
Johnson-noise thermometry has proven useful to study 
carbon nanotube a 11 ' 12 and more recently graphen o 4 ' 13 . 
In the high-temperature regime T p h > Tbg, w ^ find an 
energy relaxation rate, Jsc = A (T 3 — T^ h ), as pre- 
dicted by the supercollision mechanism- , with a prefac- 
tor A related to the amount of disorder and the car- 
rier density. Uniquely in graphene, there is a regime 
where disorder-assisted SCs dominate over the conven- 
tional electron-phonon collisions 7 "— above Tbg- In the 
low-temperature regime, for T p h < Tbg, we regain the 
dependence J oc (T 4 — T 4 ^) 4 * 7 - ' 10 ! 14 , which is the sig- 
nature of standard electron-phonon interaction in 2D 
graphene. 

Figure [5] shows the increase of electron temperature 
as a function of dissipated Joule power in a monolayer 
graphene device on a hexagonal boron nitride (h-BN) 
substrate. The average electron temperature T e is ex- 
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FIG. 1: Tunability of the Bloch-Griineisen temperature and noise thermometry setup, a, Electron-phonon inter- 
actions scatter carriers from one point on the Fermi surface (red circle) to another, within the boundary of the available phonon 
space (blue circle). In low-temperature regime (T p h < Tbg), (/max is smaller than which represents a full backscattering 
of electrons, b, The Fermi surface shrinks as the carrier density decreases, resulting in a smaller value of Tbg- Here, when 
T p h = Tbg, W just equals 2/cf- c, In the vicinity of the charge neutrality, one enters the high-temperature regime, where 
T p h > Tbg- Here, only phonons with q < 2k f can scatter off the electrons in the ordinary collisions (green arrows), whereas 
the entire thermal distribution of phonons is allowed for disorder-assisted supercollisions (purple arrow), d, Sketch of our 
noise thermometry setup. Measurements are carried out in liquid helium immersion: A Joule power is supplied to the sample, 
creating a hot carrier population. The hot electrons (diffuse red) induce a thermal noise (current generator Si) and are cooled 
by the graphene lattice (small balls). The current fluctuations are amplified using a cryogenic low noise amplifier (LNA). 



tracted from the shot noise Si, measured in the MHz to 
GHz band to overcome the environmental 1// contribu- 
tion (see Methods). With the application of the Joule 
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FIG. 2: Electron temperature in graphene as a func- 
tion of supplied Joule power. Electron temperature as a 
function of Joule power per unit area, P = V 2 /(RWL), in a 
sample of width W = 2.8 fim and length L = 2.2 /im. Data is 
shown for selected gate voltages: V s ranges from high carrier 
density (-55 V) to charge neutrality (V g = 12 V). For compar- 
ison, the lattice temperature (dashed line) is estimated from 
T ph « (P/Eif) 1 ' 4 with the lattice-substrate coupling constant 
Ex — 10 Wm~ 2 K -4 . Inset: Resistance as a function of gate 
voltage at 4.2 K. 



power P to the electron system, T e rises well above the 
substrate temperature To = 4.2 K: T e reaches 300-700 K 
when ~1 mW of Joule power is generated in our ^5 /im 2 
sized sample. The large increase of T e is, partly, due to 
the much smaller carrier density in graphene compared 
to a conventional metal. The Joule power is dissipated 
in a relatively small number of carriers, resulting in the 
corresponding increase of T e . Accordingly, we find T e 
to be at its largest value at the charge neutrality point 
(CNP), located at V g = 12 V (inset of Fig. [2]). Moving 
away from the CNP (i.e. increasing the carrier density 
n s ) T e decreases at constant power P. 

The large thermal decoupling of electrons and phonons 
(T e T p h > To) reflects the weak electron-phonon 
interaction in graphene. The steady-state value of T e 
is determined by the balance between the Joule heat- 
ing P and the cooling powers J at play—. The cool- 
ing of electrons occurs via either heat transfer to the 
phonons or heat diffusion to the metallic leads. At suf- 
ficiently high bias the diffusive contribution can be ne- 
glected, and the electron-lattice cooling can be directly 
investigated 4 -. Similarly, T p h reaches a steady-state value 
when the same power from the electrons is transferred 
to the substrate with a typical black-body radiation law 
P = Sif(T 4 ^ — TJ 4 ). The phonon-substrate coupling 
constant 16 , Ek — 10 Wm~ 2 K~ 4 , is 3-4 orders of mag- 
nitude larger than the electron-phonon coupling con- 
stant in graphene 4 -, and we achieve phonon temperatures 
T ph ~ 4-65 K with P = 0-0.2 mW/^m 2 . 

In order to test the supercollision mechanism, we plot 
in Fig. [3] the electron temperature in the representation 
Tg/P as function of P. SC theory predicts the energy 
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FIG. 3: Super-collision cubic law. Electron temperatures 
arc plotted as T'i / P for a range of selected gate voltages. Tbg 
is tuned by the electrostatic gate potential, and displayed by 
the dashed-dotted line. For temperatures T p h > Tbg, the 
SC cooling dominates, visible from the arising plateaus in the 
Tg/P representation. 
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while modeling disorder by short-range scatterers. Here g 
is the electron-phonon coupling, v is the density of states 
per spin/valley flavor, and / is the mean free path. In 
this model kpl is a constant inversely proportional to the 
strength and concentration of impurities. Other types 
of disorder like ripples 1 * or to some extend long range 
impurities can give rise to a similar T 3 dependence with 
however different expressions for the coupling constant. 
Note that, with T e 3> T p h, the T 3 /P plateau corresponds 
to the inverse of the coupling constant A. 

The T e 3 /P plots (Fig. [3]) demonstrate three different 
cooling regimes; the dip at P representing the heat 
conduction to the leads, the asymptotic P oc T 3 behavior 
(the T 3 / P plateau) pronounced in the high-temperature 
regime (T p h > Tbg), and the low-temperature P oc T 4 
behavior— at T p h < Tbg which translates into a decreas- 
ing function T 3 /P ~ P -1 / 4 . In order to highlight the 
border between the low and high temperature regimes, 
the gate-tuned Tbg is visualized by a dashed-dotted line 
which connects the values of Tbg at each V g , while T p h 
is denoted in the upper x-axis of the Fig. [3l We use 
the relation Tbg = {2v s /vf)Ef /Jcb ~ 54^^ K where 
the carrier density n s = \J[C g {V g — VcNp)/e] 2 + «o ( m 
units of 10 12 cm~ 2 ) is deduced from gate voltage V g using 
the gate capacitance C g = 35 aF/fim 2 . C g is estimated 
from device geometry. We have introduced a residual 



density no~4xl0 11 cm 2 to account for the minimum 
conductivity due to electron-hole puddles at the CNP. 

Our observation of the T 3 /P plateaus, which are most 
pronounced near the charge neutrality, strongly sup- 
ports the SC cooling mechanism in the high-temperature 
regime. With increasing n s (resp. Tbg oc y 7 ^!), plateaus 
start to develop at larger P (resp. Tp^), up to a point 
where the low-temperature regime (Tph < Tbg) domi- 
nates over the entire power range (for V g ^$ — 32 V). 

We now focus on the high-temperature regime and 
compare in greater detail the SC theory to our data. 
The short-range disorder-assisted SCs, shown in Eq. ((T|), 
scale linearly with carrier density, n s ~ v 2 (fi). In Fig. 21 
the SC coupling constant A, extracted from the T 3 /P 
plateau, is presented for — 20 V < V g < 12 V (red dots). 
Indeed, we find a linear dependence A ot n s , in a good 
agreement with the SC theory. Evaluating Eq. (JXJ) , we 

obtain A 7.5 x 1CT 4 ^f- [W/(m 2 K 3 )], where D is 
the deformation potential in the unit of eV, and n s in 
the unit of 10 12 cm -2 . Taking kpl = o-J^ ~ 3.5, from 
conductivity data, the slope in Fig. @] indicates a defor- 
mation potential D ss 70 eV, larger than the reported 
values D ~ 10 - 30 eV^iLiS. We are led to the conclu- 
sion that although Eq. (JT]) gives the order of magnitude, 
correct carrier density, and temperature dependence of 
the cooling power, it also leads to an underestimate of 
the coupling constant. This is possibly an indication that 
the short-range scattering hypothesis is too restrictive. 
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FIG. 4: SC Coupling constant A as a function of car- 
rier density. The coupling constant A increases linearly with 
the carrier density and is dependent on the level of disorder 
in the sample. The latter is changed by a heat treatment: 
red dots correspond to data acquired before and blue squares 
to measurements after the annealing. Right inset: Coupling 
constant A as a function gate voltage, deduced from the T^/P 
plateaus for steps of 2 V in the range of -20 V < V g < 12 V. 
Left inset: Resistance versus (V g — Vcnp) measured at 4.2 K 
before (red) and after (blue) the heat treatment. 
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In order to check the disorder dependence in A, we 
have increased the impurity content of the sample by 
a heat treatment; this has shifted the CNP to 115 V, 
and decreased the mobility accordingly (see left inset of 
Fig. [4]). Consistently we observe a ~ 30% increase of 
the the coupling constant (blue squares in Fig. SJ, which 
confirms the effect of impurity concentration as predicted 
in Eq. ©. 

Turning to the discussion of the results, the enhance- 
ment factor for the disorder-assisted SC cooling over the 
conventional cooling pathways, J oc (T e — T p h), is given 
as m Ref . [0 

Jsc ^ 077 T 2 +T e T ph +T 2 h 
J k F l T| G 

At the CNP, we have T e ss 400 K and T ph ~ T BG ~ 40 K 
for P = 0.025 mW/^m 2 (Figs. [2] and The enhance- 
ment factor Jsc/Jq is as large as 80 times, and explains 
the immediate observation of the SC cubic law right 
above the Tbg- Note that with further increase of the 
P (T e and T p h), the enhancement factor becomes even 
larger. 

While the SCs cubic behavior and the linear law 
for normal collisions are valid in the degenerate limit, 
ksT e < Ep, the normal clcctron-phonon collisions are 
predicted to give J oc T 4 (T e - T ph ) for k B T e > Ep&lS.. 
In our experiment, the residual carrier density at the 
CNP limits Ep above 65 meV, and we could not achieve 
the regime ksT e 3> Ep. Higher mobility samples are 
needed to access this non-degenerate regime as well as 
electron cooling by optical phonons which are both of 
great interest. 

In conclusion, we have investigated impurity mediated 
electron-phonon interaction in diffusive graphene by mea- 
suring the energy-loss of hot-electrons both below and 
above the Bloch-Gruneisen temperature for phonons. We 
observe the T e 3 dependence on the electronic tempera- 
ture predicted for supercollisions at high temperature 
while recovering the T 4 dependence for ordinary colli- 
sions at low temperature. The Bloch-Gruneisen crossover 
temperature agrees with the estimated phonon temper- 
ature and the coupling constant is consistent with the 
carrier density and disorder dependence predicted for 
short-range impurity scattering. Beside its implication 
for electron-phonon physics, our work is of direct rele- 
vance for the performance of graphene bolometers and 
photo-detectors^—. 

During the writing of this manuscript, we have become 



aware of the preprint 2 ^ dealing with the SC cooling in 
graphene. In this experiment the photocurrent generated 
at a graphene p-n junction is well described by the SC 
cubic law. 



I. METHODS 

The experiments have been performed on exfoliated 
monolayer graphene on h-BN/Si02/Si substrates. The 
h-BN platelet was exfoliated from a high quality h-BN 
powder (St. Gobain "Tres BN" .) The graphene flake was 
subsequently placed on top of the h-BN using a transfer 
technique described in Ref. [U The heavily p-doped Si 
was used as a back gate and the thickness of the insulat- 
ing layer was ~1 /im. The samples, produced by means 
of e-beam lithography and dry etching, were embedded 
in a coplanar wave-guide adapted for GHz frequencies 25 . 

We applied a Joule power P = V 2 /R to create the hot 
carrier population and study the electronic cooling as 
shown in Figs. [2] and [3] The sample's length and width 
were L = 2.2 ^m and W — 2.8 fim, respectively. The 
sample resistance R ~ 1 /n s was calculated from the volt- 
age drop across the bias resistance Rbias = 4.7 kfi (see 
Fig. HJl) . The carrier density n s was controlled electro- 
statically by means of the back-gate voltage V g . 

All data was measured for samples in a direct contact 
with liquid helium to ensure a cold phonon bath. The 
noise signal resulting from the hot carriers was first am- 
plified by a cryogenic low noise amplifier at T = 4.2 K 
followed by two additional amplifiers at room tempera- 
ture (overall gain of 82 dB). A 460 Ohms Al/AlOx/Al 
tunnel junction was used to calibrate the shot noise. Our 
broad-band setup, with a bandwidth of ~1 GHz, allows 
us to quantitatively separate the white shot noise con- 
tribution Si from 1// noise by fitting the spectra with 
a Si + C / f law. We find C cx I 2 in accordance with 
the Hoodge law 2 ^. This procedure is especially impor- 
tant for small samples and high currents where the 1/f 
contribution can be large. 
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